Quasi-bound states of quantum dots in single and bilayer graphene 
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Dirac fermions interacting with a cylindrically symmetric quantum dot potential created in single 
and bilayer graphene are not confined but form quasi-bound states. The broadening of these quasi- 
bound states (i. e. the inverse of their lifetimes) decreases (increases) with the orbital momentum 
of the electron in the case of graphene (bilayer). Quasi-bound states with energy below (above) 
the barrier height are dominantly electron(hole)-like. A remarkable decrease of the energy level 
broadening is predicted for electron energies close to the barrier height, which are a consequence of 
the total internal reflection of the electronic wave at the dot edge. 



PACS numbers: 73.63.Kv, 73.43.Cd, 81.05.Uw 



I. INTRODUCTION 

Quantum dots, or also called "artificial atoms", [l[ are 
one of the most intensely studied subjects in present day 
condensed matter physics. Recently, the realization of 
stable single layer and bilayer carbon crystals (graphene) 
has aroused considerable interest in the study of their 
electronic properties 0, [H, 0, H, @ . These new systems 
exhibit special excitations which are described by the 
analogs of the relativistic Dirac equations 0]. One of 
the most impressive effects is the so called Klein para- 
dox according to which electrons can cross large electric 
barriers with unity probability [H, § . As a consequence 
the control of the electron behavior by means of electri- 
cal potentials and structures becomes a very challenging 
task and the creation of quantum dots in such materials 
is not obvious. 

It follows from the Klein effect straightforwardly that 
the electron will escape from any potential minimum, and 
that there are no bound states in an electrically defined 
quantum dot (except in gated bilayer[10]). However, the 
polar diagrams presented in Ref. [8j indicate that in sin- 
gle and bilayer graphene the electron penetration into a 
potential barrier is strongly reduced if this electron prop- 
agates at some angles with respect to the barrier. From 
this observation we may expect that there might be long 
living quasi-bound states in such quantum dots for spe- 
cific orbital momenta of the electrons. These states may 
be probed in experiments, say by tunneling currents di- 
rected perpendicular to the dot using e. g. STM, or in 
the near field infrared absorption, as narrow peaks in the 
local density of states. 

In this work we impose a circularly symmetric quan- 
tum dot potential as created e. g. by the split gate tech- 
nique in single layer and bilayer graphene and discuss the 
conditions under which quasi-bounded states can appear. 

Such a problem was recently discussed in Ref. [ll| 



where the semiclassical approach was applied to the prob- 
lem of electron motion in a parabolic quantum dot and 
the imaginary part of the energy eigenvalues (i. e. the life- 
time of those states) was calculated. The possibility to 
confine the electrons by an external potential in a small 
region of a graphene strip was also discussed in Re f. 1131 . 
Our approach is essentially different from Ref. [ll| in 
two ways: 1) wc define the quasi-bound states through 
the averaged local density of states, and 2) we study a 
quantum dot with a step profile instead of a parabolic dot 
where the potential tends to infinity at large distances. 
We found that the width of quasi-bound states (i. e. the 
inverse of their lifetime): 1) has the opposite dependence 
on the angular quantum number for a single and a bi- 
layer graphene and 2) that it becomes extremely small 
for energies near the potential barrier height. 

The paper is organized as follows. In section II wc 
formulate the problem of a quantum dot in graphene. 
Its solution is given in section III, and in section IV the 
averaged local density of states is considered. The results 
for a single graphene layer are presented in section V. In 
section VI the same quantum dot problem is formulated 
for a bilayer of graphene, and its solution and results are 
given in section VII. Our conclusions are formulated in 
section VIII. 



II. QUANTUM DOT IN GRAPHENE 

From the point of view of its electronic properties, 
graphene can be considered as a two-dimensional zero- 
gap semiconductor with its low-energy quasiparticles 
(electrons and holes) described by the Dirac-like Hamil- 
tonian M: 



Ho 



(1) 
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where vp f=a 10 6 ms _1 is the Fermi velocity and cr = 
(<r x ,a y ) are the Pauli matrices. We assume that the to- 
tal Hamiltonian consists of the above Hamiltonian of free 
particles in addition to the cylindrically symmetric elec- 
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trie confinement potential 



V(r) = VQ(a-r) 



0, < r < a; 
V, a r < oo. 



(2) 



Such potential can be created by means of a patterned 
gate electrode with the edge smearing much less than 
the characteristic Fermi wavelength of the electrons and 
in turn much larger than the graphene lattice constant. 

As the electric potential can not confine the electrons 
in a finite region of the graphene plane there are no bound 
states. Consequently, the electron will have a continuous 
spectrum, and electron states in the quantum dot have to 
be described as decaying quasi-bound states. Neverthe- 
less we shall consider this problem as a stationary one by 
artificially confining the electron within a large sample. 
We shall solve the stationary equation 



{H-E}^(r) = 0, H = H a + V(r) 



(3) 



in a finite graphene circle of radius R. In this case the 
energy spectrum consists of discrete levels separated by 
intervals which go to zero as when R tends to infin- 
ity. The presence of quasi-bound states in the quantum 
dot can show up as peaks in the averaged electron den- 
sity in the dot. They can be revealed experimentally as 
peaks in the tunneling current through this dot in e. g. a 
STM experiment. 



III. SOLUTION OF EIGENVALUE PROBLEM 

For the sake of convenience we introduce dimension- 
less variables, measuring the distance in the units of the 
quantum dot radius a, and energies in Hvf/ci units. For 
instance, for a dot with radius a = 0.1 /zm the above en- 
ergy unit is of 6 meV. It enables us to write the Hamil- 
tonian as 

/ VQ(l-r) -i(d/dx-id/dy)\ (] 
\-i(d/dx + id/dy) ve(l-r) )' (> 

Lets write the two component wave function as 



*(r) = 



A(r) 
B(r)J ' 



(5) 



and taking into account the cylindrical symmetry of the 
problem, where we use 

^--i± = e-^(-- 1 -—) (6) 
dx dy \dr r dip J ' 

the eigenfunction problem ([3]) can be written as two cou- 
pled differential equations for the wave function compo- 
nents: 

V-^ A = ^(*-^) B > (?a) 
(V-^ = te* + A (7b) 



These equations can be further simplified by using the 
circle symmetry of our problem which allows us to assume 
the following angular dependence of the wave function 
components: 



A(r) 
B(r) 



a{r) 
ie ltp b(r) 



(8) 



where the integer m stands for the eigenstate angular 
momentum. This assumption converts Eqs. (JT)) into the 
following set of coupled ordinary radial differential equa- 
tions: 



(V-E)a 
(V-E)b 



d m + 1 

dr r 
d m 
dr 



b, 



(9b) 



These two equations have to be solved in the inner (r < 1) 
and outer (1 < r < R) region of the dot ensuring the 
continuity of both wave function components (a and b) 
at the quantum dot edge r = 1. Moreover, the proper 
boundary condition has to be satisfied at the sample edge 
(r = R). Although the exact boundary condition de- 
pends on which sublattice atoms are on the sample edge, 
we restrict our consideration to the simple equation 



a(R) = 0, 



(10) 



as the average local density of states which we are looking 
for is not sensitive to the microscopic details of the sample 
edge. 

Now inserting b from Eq. (|9bj) into Eq. ([9a| we arrive 
at the second order ordinary differential equation 



dr 2 



1 d 

r dr 



-(V-Efa, (11) 



which actually coincides with the Bessel function equa- 
tion. The other wave function component can be easily 
obtained from Eq. (f9b|) . Thus, the solution inside the dot 
(where V(r) = 0) in the case of positive energy E > 
reads 



a = FJ m (Er), 
b = FJ rn+1 {Er). 



(12a) 
(12b) 



Note we did not include the Bessel function of second 
order Y m (Er) into our solution, as it is singular at r = 0. 
Outside the dot (1 ^ r < R) the solution is 

a = PJ m (nr) + QY m (Kr), (13a) 
b = T{PJm+i(Kr)+QY m+1 ( K r)}, (13b) 

where 



k = \E-V\, 



(14) 



and the sign in the right hand side of the b expression 
coincides with the sign of (E — V). 
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Satisfying the boundary conditions we obtain the fol- 
lowing set of algebraic equations for the coefficients (F, 
P, Q): 

FJ m {E) = PJ m {K) + QY m ( K ), (15a) 
FJ rn+1 (E) = t {PJ m +i(K) + QY m+1 (K)} , (15b) 
PJ m {nR) + QY m {nR) = 0. (15c) 

As we are interested in the limiting case R — ► oo, we 
replace the Bessel functions in Eq. (|15cp by their asymp- 
totic expressions, which results into 



P cos ( kR 



nm n\ „ / nm n\ 

— --)+Q S u.(kR- — - i )=0. 

(16) 

Up to a normalization factor the solution of this equation 
can be chosen as 



P 



sin (kR 



nm 



= — cos ( kR 



TV 

4 

nm 7T 

~Y~ ~ 4 



(17a) 
(17b) 



IV. LOCAL DENSITY OF STATES 



Inserting the obtained solution for the P and Q coeffi- 
cients into Eqs. (|15h .b) we arrive at the equation for the 
eigenvalues. It follows from Eqs. (flT|) that these eigen- 
values are separated by 



AE = An = —. 

R 



(18) 



These dense discrete energy levels are a consequence of 
the finite size of our sample and are not of interest to us. 
Therefore we choose another procedure and solve the two 
equations (fTSh .b) together with the condition 



P l 



Q =1 



(19) 



which is clearly satisfied as it follows from Eq. (p~7|) . This 
procedure enables us to obtain the local density of states 
from which we can derive the quasi-bound states. The 
procedure is as follows. 

First, we calculate the normalization factor N of the 
wave function. In the limiting case R — > oo it can be 
done just using the asymptotic wave function expression, 
namely, 



1 = 2- 2nN z 



rdr 

V 7TKV 



jp COS (f 



■nm 7r\ 1 2 4:N 2 R . 9 9N 



(20) 



The additional factor 2 appears because both wave func- 
tion components have to be taken into account. Using 
Eq. we find 



N 



(21) 



Next, we point out that any physical property of the 
quantum dot, say like the tunneling current through the 
dot, or the absorption of the infrared radiation in near 
field spectroscopy, can be expressed as a summation of 
some matrix elements over the above dense quantum 
states. The matrix elements are integrals over the quan- 
tum dot area, namely, 



M = 2ttN 2 F 2 



rdrf(r){j 2 n (Er) + Jl +1 (Er)} 



(22) 

where the function f(r) characterizes the interaction of 
the quantum dot with the measuring probe. Thus, re- 
placing the summation over the discrete levels by an in- 
tegration, in accordance with Eq. (TT5)) . 



^2 M w — J dEM - 2i? J dEN 2 F 2 = J dE 



kF 1 



we see that within the accuracy of the experimental form- 
factor (the integral in the right hand side of Eq. (|2"2")l ) the 
quantity 



\E-V\F 2 



(24) 



acts as a local density of states in the quantum dot area. 



V. RESULTS FOR QUANTUM DOT IN 
GRAPHENE 

We solved Eqs. (p~5b .b fT^|) numerically from which we 
obtained the coefficients, and the local density of states 
([24]) . A typical example for the two components of the 
wave function together with the confinement potential 
profile is shown in Fig. [T] 




FIG. 1: (Color on line) Profile of the confinement potential 
(dot-dash curve), and the two wave function components: a - 
solid curve, b - dashed curve. Barrier height V = 10, energy 
E = 6, and orbital momentum m = 2. 



As the energy is lower than the potential height we 
see that the two wave function components have a differ- 
ent phase indicating the electronic type character of the 
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wave function inside the dot, and the hole type character 
outside it. The large value of the wave function com- 
ponents inside the dot show that this eigenfunction is a 
quasi-bound state. 

Typical local density of states are shown in Fig. [2j It 
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FIG. 3: (Color on line) Quasi-bound states with orbital mo- 
mentum m = for a quantum dot in graphene. The energy 
of these states are given by the black curves and its width 
(i. e. the inverse of the life time) by the colored region. The 
straight slanted line corresponds to E = V. 



FIG. 2: (Color on line) Local density of states as function of 
the energy in case of barrier height V = 12 for the orbital 
momenta m = (solid curve), m = 1 (dotted curve), and 
m = 2 (dashed curve). 



exhibits peaks which can be associated with the quasi- 
bound states of the dot. The three curves correspond to 
the following orbital momenta of electron m — 0, 1, 2. We 
observe the general tendency that the larger the orbital 
momentum the narrower the peaks. Noticeable very nar- 
row peak when the energy is close to the barrier height 
(see the curve for m — 0). This tendency is even better 
seen in Figs.[3]and|4]where the positions and broadenings 
of the peaks are shown. 

We fitted peaks in the density of states by Lorentzian 
functions a n j n /{(E — E n ) 2 + 7„} defining three param- 
eters for any of them: the position E n , its broadening 
7„, and the amplitude a n . Graphically these parameters 
are shown in Figs. [3] and [4] for two orbital momentum m 
values as function of the barrier height V. The posi- 
tions E n of the quasi-bound levels are shown by the solid 
curves while the shaded areas between two i5„±7„ curves 
indicate the broadening of them. 

As expected, in the case m — the levels are rather 
broad. Actually, they can hardly be identified as quasi- 
bound levels and they rather correspond to weak oscil- 
lations in the local density of states of the continuous 
spectrum (see the solid curve in Fig. [J). 

In the case of to = 2 we see (Fig. 2]) a quite different 
picture. The levels are narrow and indicate the presence 
of long living quasi-bound states. It is interesting to see 
that the quasi-bound states are seen above as well as 
below the barrier, the latter is indicated by the slant solid 
line. Actually, this is the consequence of the equivalence 
of the Dirac electrons and holes in the barrier region. 

As was already mentioned in Figs. [3] and 2] we see one 
more important peculiarity of the local density of states 
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FIG. 4: (Color on line) The same as Fig.[3]but now for m — 2. 



of a quantum dot in graphene. This is the extremely nar- 
row states in the vicinity of the top of the barrier. This is 
not an accidental phenomena, but the consequence of the 
important fact that it is rather difficult for Dirac elec- 
trons to penetrate the barrier when its energy is close 
to the barrier height. This property follows straightfor- 
wardly from the electron penetration through the barrier 
problem solved in [1, [II , although they payed no atten- 
tion to this limit case. The matter is that the angle <p 
(with respect to the perpendicular to the barrier, see the 
inset in Fig. [5]) of the incident electron and the angle ip 
of the refracted electron has to satisfy the equation 

Esm<p = (V-E)smi/;, (25) 

which is the equivalent of Snell's law in optics [13j |. In 
the case when the electron energy is close to the barrier 
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height 



E = V 



|A| « V, 



(26) 



the electron wave goes from the material with large re- 
fraction index into the material with small refraction in- 
dex. In this case the well known phenomena of total 
internal reflection takes place. It means that there is a 
critical incident angle 



= |A|/V, 



(27) 



such that electrons with larger incident angles (\(p\ > cpo) 
are totally reflected from the barrier (see the inset in 
Fig. El where the angles at which the electron penetrates 
the barrier are shown by the shadowed sector). The elec- 
tron current (which can be named as tunneling proba- 
bility) in the barrier perpendicular to the barrier edge 
direction can be expressed as follows: 



W 



(28) 



and is shown in Fig. [5] 




<p/cp 



FIG. 5: (Color on line) Tunneling probability dependence on 
the incident angle of the electron beam in the case of small 
deviation of the electron energy from the barrier height. Inset: 
definition of the critical angle <^o- 



So, when the electron energy is close to the barrier 
height there is only a very small region of the incident 
angles at which the electron can penetrate the barrier. It 
means that for these electrons the electrical barrier acts 
almost as a confining potential. This explains why in this 
energy region very narrow quasi-bound levels appear even 
for zero orbital momentum. 



VI. QUANTUM DOT IN BILAYER 

It is instructive to consider the same quantum dot 
problem for bilayer graphene. There are similarities and 
differences between the two graphene systems [l4| • Both 
systems have gapless electron and hole spectra, but are 
described by different Dirac type Hamiltonians. 



Following Ref. Q we use the two component wave func- 
tion approximation for the bilayer which is described by 
the following Hamiltonian: 



_ 1 / {d/dx-id/dyf 

~ 2 Ud/dx + id/dy) 2 



(29) 



where the distances are measured as before in units of the 
quantum dot radius a, and the energies are measured 
in h 2 /m*a 2 units. For example, for a dot with radius 
a = 0.1 /im, and the effective mass [3] m* — 0.038m e 
the above energy unit is 0.2 meV. 

Assuming expressions for the wave function analogous 
to those for graphene ([5]) we obtain the following set of 
equations for the wave function components: 



(E-V)A 
(E - V)B 



.d_ 

dy 

. d 



2 \dx dy 



1 (d_ 

2 \dx~ 
1 ( d 



A. 



(30a) 
(30b) 



The components A and B have to satisfy the continuity 
condition together with continuity of their first radial 
derivatives at the dot edge. 

In the dot and as well outside it these two equations 
can be transformed into a single equation for any wave 
function component: 



~V 4 -{E-Vf\A 



-V 2 + (E-V) 



-V 2 - (E- V) \A = 



(31) 



Using the axial symmetry and assuming the following 
angular dependence: 



A(r) 



a imip 



a(r), 



(32) 



this component has to satisfy any of the following radial 
equations: 



dr 2 



1 d 

r dr 



,2 1 



±(E-V)-- T 



0. 



(33) 



It is evident that they are the equations for the Bessel 
and modified Bessel functions. Inside the quantum dot 
(r ^ 1) the solution is: 



FJ m {kr)+GI m (kr), 



(34) 



where k = y2E. Two other functions (Y m and K m ) 
are not included because of their singularity at the origin 
r = 0. Similarly, using Eq. (|30b[) we obtain for the other 
component of the wave function 



B(r) = e i(m+2 H(r), 

b = FJ rn+2 (kr) + GI m+2 {kr). 



(35a) 
(35b) 
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Outside the quantum dot the wave function compo- 
nents are given by 



A 
D 



PJ m {nr) + QY m {nr) + SK m ( 
T{PJm+2(Kr) + QY m+2 {nr) 
+ SK m+2 {nr)}, 



Kr) 



(36a) 
(36b) 



where k = y/2\E — V\, and the sign of the B component 
coincides with the sign of the expression (E — V). We do 
not included the functions /„ and I m +2 into the above 
expressions, as these functions are responsible only for 
satisfying the boundary condition at the sample edge r = 
i?, and are therefore not relevant in the limit R — > oo. 

In analogy to our approach for graphene we equate 
both wave function components and their derivatives at 
the quantum dot edge (r = 1), and together with Eq. (fl9|) 
we obtain a set of five algebraic equations for the five 
parameters F, G, P, Q and S. 

We point out that the functions I m (kr) and K m {nr) 
together with the analogs for orbital momentum m + 2 
are essential only in the region close to the dot edge, 
where they ensure the continuity of the derivatives of the 
wave function components. While the main contribu- 
tion to the local density of states is determined by the 
functions J m (kr) and J m +2 (&?*)■ Following the procedure 
presented in Sec. IIVI we obtain the following expression 
for the quantum dot local density of states in a bilayer: 



P{E) 



1 



-F* 



(37) 



Note that it differs from the analogous expression for the 
case of graphene by the factor depending on the electron 
energy in the barrier, what is caused by the different 
dispersion law in both materials. 



VII. RESULTS FOR THE DOT IN BILAYER 

Solving numerically the set equations for the coeffi- 
cients, introduced in the previous section, we calculated 
the wave functions and the local density of states. An ex- 
ample of the wave function is shown in Fig. [5] Comparing 
it with the wave function for the dot in graphene we 
notice several differences. The derivative of the wave 
function in a bilayer is continuous at the dot edge, while 
in graphene the wave function components exhibit kinks 
there. That is the reason why in a bilayer there is an in- 
termediate region of exponential behavior (see the shad- 
owed rectangle in Fig.[S|) where the electron type function 
changes itself into the hole type one. 

The local density of states of the quantum dot states 
is shown in Figs. [7] and [8] in the case of two orbital mo- 
menta (m = and m=2). These pictures also differ 
essentially from those for graphene shown in Figs. [3] and 
|U First, we see that the quasi-bound states in a bilayer 
are much narrower as compared with these in graphene. 
The main reason of this difference is as follows. Although 
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FIG. 6: (Color on line) Wave function components in bilayer: 
V = 20, E — 10, m — 0: solid curve - component a, dashed 
curve - component b, the dash-dotted curve is the potential 
profile of the dot. 




FIG. 7: (Color on line) Quasi-bound states for quantum dot 
in bilayer with orbital momentum m — 0. The energy of these 
states are given by the black curves and its width (i. e. the 
inverse of the life time) by the shadowed regions. The straight 
slanted line corresponds to E = V. 



single layer graphene and bilayer graphene are both gap- 
less materials, the physical nature of their bands is quite 
different. In graphene the electron and hole parts of the 
bands are the natural prolongation of each other. Elec- 
trons and holes are like different expressions of the same 
Dirac quasi-particle. While in bilayer graphene these two 
contiguous electron and hole bands are much more inde- 
pendent from one another, and reminds one to the acci- 
dental touch of their extremum points. 

This can be more clearly demonstrated by rewriting 
the component equations (|30p in the case of perpendicu- 
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FIG. 8: (Color on line) The same as Fig.[7]but now for m = 2. 



lar incidence of the electron to the barrier: 

{ E-V)A=\£ B , (38a) 

(E-V)B=l^A. (38b) 

Adding and subtracting these equations we obtain the 
uncoupled equation set for the hole A + B and electron 
A — B type wave functions. Consequently, in the case 
of perpendicular incidence the electron in a bilayer can 
be confined in the dot by the electrical potential, and 
thus, we obtain stationary states. In the case of slanted 
incidence, when the electron momentum component k\\ 
along the barrier is not zero the equations can no longer 
be decoupled. In this case this longitudinal momentum 
(or the orbital momentum of electron) serves as a cou- 
pling constant between the electron and the hole. Con- 
sequently, the larger the angular momentum, the larger 
the probability for the electron to convert itself into a 
hole, or the more dominant is the Klein effect, and thus 
the more smeared is the quasi-bound level. This is ac- 
tually seen in Figs. [7] and [5] which show that the quasi- 
bound states for m = are narrower as compared to the 
to — 2 states. By the way, this difficult penetration of 
the electron into the barrier is indicated by the interme- 
diate exponential region in the wave function coordinate 
dependence as shown by the shadowed area in Fig. [6] 

This simple physical picture explains one more inter- 
esting property of the above local dot density of states. 
In Figs. [3] and 2] we see that in graphene there is some 



symmetry between the quasi-bound states below the top 
of the barrier (blue solid slanted line) and above it. While 
in bilayer graphene as seen in Figs. [7] and [8] such symme- 
try is absent and the states above and below the top are 
quite different. The states above the top of the barrier 
are much more smeared. This is caused by the fact that 
above the barrier the nature of the wave function is the 
same in both regions (in the dot and outside it), and 
there is no need for the electron to transform itself into 
a hole, and consequently, the probability to escape the 
dot is larger than in the opposite case when the energy 
is smaller than the top of the barrier. 

VIII. SUMMARY AND CONCLUSIONS 

Using the two wave function component approxima- 
tion we calculated the local density of states in an elec- 
trically defined circle symmetric quantum dot in single 
layer and bilayer graphene. It was shown that in bilayer 
rather narrow quasi-bound states appear when the en- 
ergy is smaller than the barrier height. The broadening 
of the states in bilayer graphene increases as the orbital 
momentum becomes larger which is opposite to the case 
of graphene. 

In contrast in graphene narrow quasi-bound states are 
predicted with increasing orbital momentum. This differ- 
ent behavior of the quasi-bound states in graphene and 
bilayer is explained by the different physical nature of the 
touching electron and hole energy bands. 

Weakly broadened quasi-bound states are predicted in 
both graphene and bilayer graphene in the region where 
the electron energy is close to the top of the barrier. This 
phenomena can be understood from an analog of the op- 
tical effect of total internal reflection which an electron 
wave suffers in the above mentioned region of energies. 

We also notice the different symmetry of the states 
above and below the top of the barrier in single layer 
graphene and bilayer graphene, caused by the different 
way of electron conversion into a hole in both systems. 
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